clear all
clc 
i=1;
for  k=17.5:0.01:17.8
    m=1.109*10^5;
    k3=0.42*k;
    a=29.7*10^9*0.5*exp(-1/3*(4+2*k3))+10.1*10^9*exp(-(4/3+2/3*k3+2/3*k))...
        +3.29*10^9*exp(-1/3*(4+4/3*k));
    l1=1.22*10^11*m*(10.1*10^9*0.5*exp(-2/3*(0.8+k3))+...
        3.29*10^9*exp(-2/3*(0.8+k)))/(1.22*10^11+...
        m*(10.1*10^9*0.5*exp(-2/3*(0.8+k3))+3.29*10^9*exp(-2/3*(0.8+k))));
    l2=1.22*10^11*m*a/(1.22*10^11+m*a);
    l3=10.1*10^9*exp(-(0.8+2/3*k3))+3.29*10^9*exp(-(0.8+2/3*k));
    y=fun(k3,k,m,l1,l2,l3);
    K(i)=k;
    Y(i)=y;
    i=i+1;
end

[m,n]=size(K);
i=1;
j=2;
for i=1:n
    if Y(j)>Y(i)
        j=j+1;
    else
        j=i;
       break
    end
end
K(j)
Y(i)
